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Desde que se ha hecho evidente que existe una función de 
estado p(t) que guía la evolución de las plumas de solutos 
en los flujos naturales, la naturaleza de esta función ha sido 
analizada desde diversos puntos de vista (Constaín, 2012). 
En este artículo se muestra cómo q(t) es universal para todo 
el flujo en tanto éstos cumplan la condición de “equilibrio 
dinámico”. Al mismo tiempo se muestra la validez 
simultánea de la ecuación de Elder en función de q(t), que 
relaciona el transporte dispersivo con la geomorfología y 
por lo tanto pudiendo calcularse la pendiente de un flujo 
en “estado estable”. Se aplica esta metodología a un cauce 
mediano en Colombia, verificándose que la gráfica de la 
concentración pico para el trazador es bastante congruente 
con los datos esperados si se aplican los principios aquí 
desarrollados. 


Palabras clave: estudios de calidad de agua, hidráulica, 
dispersión, geomorfología. 


Introducción 


Un hecho observado en los procesos fluviales 
es que existe una relación evidente entre las 
fuerzas que tallan los cursos naturales de agua 
y la forma y proporciones de los lechos de los 
mismos. Así, por ejemplo, en un cauce estrecho 
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la fuerza del agua tallará de manera preponde- 
rante las paredes laterales (A), tratando de ganar 
más ancho; por otro lado, en cauces anchos, esta 
fuerza tenderá a profundizar el cauce (B). De esta 
forma, los procesos de “tallaje” en un sentido u 
otro buscarán equilibrarse de modo dinámico en 
cierto punto (C) (figura 1). 
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(a) 


Figura 1. Procesos de “tallaje” en el lecho. 


Ocurre entonces un “equilibrio dinámico” 
cuando la sección transversal del cauce logra un 
valor tal que las fuerzas de “tallaje” en distintas 
direcciones no predominan la una sobre la otra 
y la geomorfología del flujo se mantiene aproxi- 
madamente en el tiempo (Leopold $: Maddock, 
1954). 

Este estado de “equilibrio dinámico” en el 
proceso de formación del cauce por el movi- 
miento continuo del agua no sólo se interpreta 
en términos de los cambios en las áreas efectivas 
de las paredes resistentes del lecho (efecto), 
sino también, e indisolublemente, mediante la 
velocidad del flujo, pendiente y rugosidad su- 
perficial (causas). El efecto final de este proceso 
de estabilización en sentidos contrarios: no sólo 
afecta la dimensión transversal sino también la 
longitudinal. Así, desde el punto de vista del 
perfil longitudinal del lecho, es que aquellos 
sitios donde haya una mayor pendiente, la 
deposición de sedimento tratará de suavizar 
el perfil y que en aquellos sitios de pendiente 
“suave”, la socavación intentará “empinar” el 
perfil, siendo Sb la pendiente (figura 2). 

El “equilibrio dinámico” en la formación 
de los cauces puede entonces verse como una 
tendencia a la igualación de las tasas de “depo- 
sición” y “socavación” de sedimentos sobre el 
lecho tanto en el perfil transversal como en el 
longitudinal. Como quiera que el agente con- 
trolador directo de estos dos procesos opuestos 
es la velocidad media del fluido, este parámetro 
no toma por supuesto valores arbitrarios sino 
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aquellos compatibles con la tendencia explícita 
al mencionado equilibrio. 

Ahora bien, para un grado dado de resisten- 
cia del lecho (dependiente de las características 
físicas de las partículas de ese lecho), esto es, 
lo que se denomina “rugosidad”, la pendiente 
como fuerza motriz del movimiento del fluido 
será un cierto valor que ajusta de forma con- 
cordante la interacción velocidad-deposición- 
socavación (Hack, 1957). 

Por lo tanto, el mencionado “equilibrio di- 
námico” se refleja sobre tasas del transporte de 
masa en el interior del sistema. Esto se analiza 
en seguida. 


Transporte de masa en el interior del 
sistema bajo la condición de “equilibrio 
dinámico” en los cauces naturales 


Considerando ahora un sistema cerrado evo- 
lucionado totalmente irreversible, (2 por el 
principio de conservación de la energía, en el 
que la energía mecánica entrante se convierte 
enteramente en energía térmica por degrada- 
ción, se tiene que: 


AU=A0Q (1) 


En este caso se puede ver que la pendiente Sh 
aporta un componente inicial de energía poten- 
cial gravitatoria, AU (Z), el cual es transformado 
en energía cinética, AK (movimiento del agua), 
el cual finalmente es transformado por entero 
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eS 


Sb 


Figura 2. Proceso de igualación del perfil. 


en calor, AQ (figura 3). 

También se ha demostrado que en los sis- 
temas hidrológicos en los que se estudian los 
cambios de la geomorfología es importante 
aplicar además el segundo principio de la 
termodinámica o principio de la entropía. En 
estas condiciones, para un sistema cerrado 
irreversible isotérmico en el que toda la energía 
potencial U se transforma en calor, Q, se puede 
mostrar entonces que la entropía es constante 
(Leopold €: Langbein, 1962). 

Entonces la entropía interna (por irreversibi- 
lidad), Si, expulsada en la frontera del sistema, 
con T' como temperatura absoluta, se puede 
escribir como: 


o (2) 


Si se parte de la definición del gas ideal para 
la nube de trazador PV = nRT, donde P es la 
presión, V es el volumen, n el número de moles 
y R =8.31 J/(K*mol), la constante universal de 


O AK 
BS 


Figura 3. Transformaciones de energía en el flujo. 


los gases, reemplazando la temperatura en la 
ecuación anterior, se tiene: 


_ (nR) AU 
PV 


ASi= (3) 

Considerando ahora la presión P = fl Af 
como fuerza, área como A = AY * AX, y volu- 
men como V = AY * AX* AZ, con Y = ancho, 


X = largo y Z = profundidad del flujo, se tiene: 


2% (nR) AU _ (nR) AU 


ASi 


Py És 
A Ay xAx 


_(MR)AU _(1R) AU 
 fxnbz u (4) 


Dividiendo cada miembro por el volumen y 
sabiendo que AV= q* At, donde q es el caudal, 


A 
A o (5) 


gq AV 


queda finalmente: 


Esto implica que la producción de entropía 
interna por unidad de caudal lleva de continuo 
ala disminución de la energía relativa disponible por 
unidad de volumen. La entropía creciente dentro 
del sistema genera una distribución cada vez 
más uniforme de energía dentro del volumen de 
O. El miembro derecho de esta ecuación puede 
ser interpretado entonces como la descripción 
de la distribución más probable de energía en el 
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sistema, la cual es constante según la interpre- 
tación estadística de Boltzmann (Leopold éz 
Langbein, 1962), en los que los distintos eventos 
energéticos en el volumen de control (2 son 
equiprobables: 


Si =kx) p, xIn(p,) = máximo (6) 
a 


PL=P2=P3=P4= =P; (7) 


El “equilibrio dinámico” (estado estable) en 
los cauces naturales va a significar entonces no 
sólo que el flujo neto de masa en las fronteras 
del sistema es aproximadamente nulo en virtud 
de la equivalencia entre socavación y deposición 
(conservación de la masa), sino también que los 
eventos individuales internos descritos median- 
te los intercambios de energía serán equivalen- 
tes, uniformes o indistinguibles en el volumen 
de control por su equiprobabilidad. 


Cálculo de las diferentes tasas de 
transporte de masa en un sistema de 
flujo en “equilibrio dinámico” 


Es interesante ahora estudiar las consecuencias 
concretas de la “equiprobabilidad” para diver- 
sos sucesos de intercambio de energía dentro del 
sistema, en este caso la traslación de diferentes 
porciones de masa dentro del sistema, sea ésta 
una masa disuelta o suspendida en un sistema 
hídrico. 

Se debe partir del hecho de que en un flujo 
natural la energía potencial estará representada 
como se indica a continuación: 


U=mxgxh (8) 


Supóngase que un cierto intercambio ele- 
mental de energía dentro de un sistema que 
por motivos de su “localidad” está asociado 
más bien a un incremento de masa que a uno 
de altura: 


AU =Amxgxh (9) 
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Si el sistema es estacionario (equilibrio di- 
námico), entonces los diferentes intercambios 
de energía derivados de las permutas de masa 
de partículas similares de soluto tendrán igual 
probabilidad. Entonces los propios intercambios 
de masa, como distribución estadística en el 
tiempo, estarán “igualados” también por este 
principio. 

Corresponde por lo tanto a una distribución 
estadística constante similar a la de la ecuación 
(5), sólo que en lugar de la variable “energía 
potencial” U se tiene la variable masa M. Inte- 
resa además establecer la distribución repartida 
además en el tiempo. 

Para analizar su distribución en el tiempo se 
define la traslación en el espacio de un elemento 
de masa Amj en un tiempo Afj. Entonces la “tasa 
de transporte de masa” es igual a la razón entre 
el diferencial de masa movilizada en un punto 
sobre el diferencial de tiempo en que ocurre este 
evento (figura 4). 

Matemáticamente, esto se define así: 


Ttm., = a, (10) 
j At, 

La tasa de transporte de masa “general” 
para el tramo considerado, definida por su 
valor representativo desde el punto de vista 
estadístico, es la esperanza matemática de 
los componentes individuales de transporte 
(Spiridonov €: Lopatkin, 1973). Por lo tanto su 
definición representativa en todo el dominio del 
volumen de control O será: 


1 Am, 
<Ttm>=* Y —— 11 
yo (11) 


Pero como los elementos discretos de masa 
transportada son equivalentes (iguales por 
equiprobables), entonces queda: 


<Ttm>=2x(¡xTtm) (12) 
) 


Finalmente: 
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Figura 4. Movilización de un elemento de masa en O. 


<Ttm>=Ttm, (13) 


O sea, la tasa de transporte de masa “gene- 
ral” es convergente con las tasas de transporte 
de masa “locales” en “estado estable” del 
cauce, consideradas ellas como indistinguibles 
(iguales). En este caso, lo general replica a lo 
particular. 


Cálculo de las diferentes tasas de 
transporte de masa en un sistema de 
flujo en “equilibrio dinámico”: caso de la 
masa disuelta 


Aunque la masa transportada en el flujo puede 
ser suspendida o disuelta, es posible en estado 
estable realizar un “balance detallado” (Prigogi- 
ne € Kondepudi, 1998) de ambos tipos de masas 
desde la óptica de la equivalencia entre “tasas 
de transporte locales” y “tasa de transporte 
general”. Ello es consecuencia del equilibrio 
dinámico en cuanto en esta condición cada 
transformación elemental es balanceada por su 
contraparte específica, en el esquema concreto 
en el que hayan sido planteadas. 
Entonces, si: 


< Ttm > (disuelta + suspendida) 


=Ttm, (disuelta + suspendida) (14) 


Se sigue que: 


< Ttm > (disuelta) = Ttm y (disuelta) (15) 
Y: 
< Ttm > (suspendida) = Ttm, (suspendida) (16) 


Por estas razones se aplicará el principio de 
equivalencia entre tasas de transporte “locales” 
y “generales” al componente disuelto, lo que 
facilita su análisis por medio de trazadores. 


Cálculo de la tasa de transporte de masa 
disuelta “general” en el sistema 


En condiciones de estado estable, en un flujo 
natural, el aporte de solutos a la corriente es una 
característica propia de la cuenca y del caudal 
al cual se incorporan estas sustancias. Si se lla- 
ma “Co” a la concentración base resultante de 
este aporte y “q” al caudal correspondiente, la 
relación típica es una curva hiperbólica (Hem, 
1985) (figura 5). 

Dicha relación tiene esta forma toda vez que 
entre mayor sea el caudal para un aporte dado 
de solutos en el lecho, se tendrá un mayor volu- 
men para disolver, esto es, menor concentración: 


q xC, = k(Cte) (17) 


Atendiendo a la parte dimensional, la “cons- 
tante” k es miligramos por segundo, o sea, una 
tasa de transporte de masa. 


Tasa de transporte de masa disuelta “local” 
medida por trazadores 


Una forma de medir las tasas “locales” de 
transporte de masa en los flujos naturales es 
por medio de trazadores, los cuales, vertidos de 
forma súbita en un punto del tramo del cauce, 
se pueden considerar como “testigos” locales 
del movimiento “individual” de la masa en el 
sistema Q (figura 6). 

Para entender de forma apropiada cómo el 
trazador puede cumplir la misión de medir la 
tasa local de transporte de masa es necesario 
estudiar algunos puntos de la dinámica de los 
mismos. 
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q(m9/s) 


Co(Mgr/1) 


Figura 5. Relación entre concentración base y caudal. 


Figura 6. El trazador como testigo “local” en el flujo. 


Velocidad media del flujo en función de los 
datos del transporte dispersivo 


En artículos anteriores, los autores han definido 
esta función de la siguiente forma: relacionando 
dos velocidades, una la de dispersión del traza- 
dor, Vdisp, de naturaleza irreversible y medida 
por su desplazamiento random walk, y la otra la 
de advección, U, que no debe confundirse con la 
energía potencial, como factor integrante. Aquí 
A y tson parámetros característicos de desplaza- 
miento y fase del movimiento Gaussiano mono- 
dimensional de la pluma de trazador (Constaín 
éz Lemos, 2011; Constaín, 2013): 
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oía 


uu u —u (18) 


Como quiera que V,,, mida la rapidez con 
que las partículas de trazador se alejan de 
modo recíproco, es una función que representa 
la reacción del sistema a la ruptura súbita del 
equilibrio químico y por lo tanto es un potencial 
termodinámico. La naturaleza especial de y(t) se 
puede caracterizar diciendo que es una función 
de estado del sistema, definida mediante la 
siguiente ecuación: 


$dp=0 (19) 


El cálculo teórico de la función q(t), como ya 
se ha explicado en otros artículos de los autores 
(Constaín $: Corredor, 2013), se hace a partir del 
tiempo medido entre los dos puntos de inflexión 
de la curva Gaussiana experimental dividida 
por el tiempo de transporte y multiplicada por 
una constante. De la ecuación (26) se puede es- 
tablecer una definición para la velocidad media 
advectiva: 


He (20) 
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Despejando E: 


pur 


E0- E 


(1) 


Debe notarse aquí que el tiempo caracte- 
rístico Gaussiano tí es diferente a la variable 
independiente t, pues la primera se refiere al 
tiempo en el cual la separación aleatoria de 
las partículas ocurre mientras esta distancia 
promedio es A, mientras que la escala temporal 
general se asocia más bien con el efecto advec- 
tivo de U. La relación entre ambos tiempos se 
puede establecer mediante una dinámica del 
tipo Poisson-Svedberg (Constaín, Peña, Mesa, 
éz Acevedo, 2014): 


0% 0,215 (22) 


La función de estado qÍt) 


Debido a que q(t) es un potencial termodiná- 
mico y en consonancia con la ecuación (19), 
su valor inicial es el mismo que el final, por lo 
tanto es una curva unimodal, con un pico en 


a”, cuya parte utilizable siempre es después de 
t=1 (figura 7). El flanco ascendente corresponde 


a los fenómenos transitorios y muy rápidos que 
ocurren cuando el soluto conservativo se acopla 
en el líquido. 


di 
yendo en la medida que el sistema tiende de 


Debido a que la velocidad V ¡sp VA disminu- 


forma irreversible hacia el equilibrio después de 
t = q, entonces q(t) es descendente en su flanco 
más externo. Debe tenerse en cuenta que hay 
valor máximo para qp(t = 0) = 2.16. 


La descripción de la pluma de soluto vista por 
diferentes observadores 


La transformación de Galileo para pasar de un 
sistema inercial a otro se incluye de manera par- 
cial en el argumento exponencial de la ecuación 
clásica de Fick, con M, la masa de soluto y A, la 
sección transversal del flujo (Fischer, 1967): 


M (uy? 
CAE 


ASES 


C(x,t)= (23) 


Sin embargo, dicha transformación es com- 
pleta sólo si el coeficiente E es una función del 
tiempo, tal como se describe en la ecuación (21). 
Ahora reemplazándola en la ecuación (23), se 
tiene con: V21tf = 1.16: 


M _ CU 
e 212 2 


2Bp 


Clx,t) = ——=— 
q, pt1.16 


(24) 


Figura 7. p(t) como función de estado unimodal. 
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La particularidad de esta ecuación (fórmula 
modificada de Fick) es que reproduce bastante 
bien las curvas experimentales de trazador, 
como se mostrará en la sección “Pantallas de 
la medición desde el equipo IDF”. Supóngase 
que se realiza un vertimiento lateral de trazador 
desde la orilla de un flujo turbulento, tal como 
se muestra en la figura 6 con la pluma de traza- 
dor evolucionando en el tramo; en este caso, el 
caudal comprendido por el “tubo de corriente” 
que limita la pluma de trazador, q, es un valor 
“local”, toda vez que, en principio, el trazador 
no ha cubierto todo el volumen del flujo real ni 
participa de la velocidad plena del flujo (0 <U) 
ni su ancho efectivo es el de todo el flujo (W, < 
W) ni tampoco que su profundidad (en la orilla) 
es la del cauce (h] < h): 


q, =W,xhxU, (25) 


Debe recordarse, además, que en un cauce 
natural turbulento vale aproximadamente una 
distribución transversal de velocidades (línea 
punteada), tal como se muestra en la figura 8. Se 
muestra también en la gráfica la anchura “local” 
de la pluma de trazador (línea roja). Se asume 
que la profundidad es el valor medio de todo el 
cauce (para simplificar, h] = h). 


Figura 8. Velocidad y ancho locales para el trazador, 


definiendo caudal local. 


+» ISSN 0187-8336 


Modelo para la evolución de la concentración 
pico del trazador en el tiempo, Cplt) 


La concentración pico del trazador, definida 
como en la ecuación (26), es muy importante, 
pues va a permitir verificar de forma global la 
validez de los principios aquí propuestos: 


M 


C,(t) = ——— 
0 q qt1.16 


Pp 


(26) 


Ahora se puede mostrar que su evolución 
temporal se puede definir más sucintamente 
mediante la siguiente definición empírica: 


C,(t) = alt) xt* (27) 


Mo.” 


El factor multiplicador “a” es una función 
del tiempo que debe ser ajustada a partir de los 
datos experimentales en diferentes instantes. 
Es conveniente partir de un modelo simple 
exponencial descendiente del tipo: 


alt)=a, e *' (28) 


Definición de la tasa de transporte de masa 


local mediante la ecuación modificada de Fick 


A partir de los conceptos desarrollados anterior- 
mente se puede ahora definir la “tasa de trans- 
porte de masa local” con base en experimentos 
de trazador. En el punto de medición para la 
ecuación (33) se tiene que la tasa local de trans- 
porte de masa es aproximadamente igual a la 
siguiente ecuación, con Cp como concentración 
pico para el trazador, considerado conservativo 
para que su disminución a lo largo del tramo sea 
debida sólo al efecto dispersivo: 


Ttm MM x xq x1.16 (29) 
SAL, ci 
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Igualación de la tasa de transporte de masa 
“general” y “local” en un flujo en “equilibrio 
dinámico” 


En correspondencia con lo establecido en el 
apartado 4, ecuación (15) para la masa disuelta 
observada de forma “general” y “local” en un 
cauce natural, se tiene que para un flujo en equi- 
librio dinámico se cumple aproximadamente 
que: 


qxC,=C, xq xq x1.16 (30) 


O sea: 


E 16 (31) 


Esta ecuación es bastante significativa, toda 
vez que q(t) es, como ya se dijo, una función de 
estado y tiene una significación termodinámica 
muy especial en tanto que refleja propiedades 
“generales” del sistema en el que evoluciona, 
aunque sea medida en zonas “locales”, así trae 
información del estado particular del proceso 
de avance de la mancha de trazador (q, y CC.) 
pero también de los procesos generales para 
el cauce en equilibrio dinámico (como q y Co). 
Por esta razón se puede considerar a la ecuación 
(31) como una definición universal para dicha 
función, es decir, válida para todo el ámbito 
del flujo en el tramo considerado. Ahora, como 
quiera que q, y C, son parámetros de la pluma 
función del tiempo, el cálculo de la ecuación 
(31) debe hacerse en un tiempo especial que 
hemos denominado “tiempo de acople”, tc. En 
este instante se puede considerar que las con- 
diciones termodinámicas predominantes en el 
estado estable hacen que la información de todo 
el proceso se refleje en un solo cálculo. Dentro 
de dicha información estará naturalmente la 
correspondiente a la pendiente en su calidad 
de causa (“driving force”) del proceso. 


Papel del “tiempo de acople”, tc, en la 
formación de la pluma de soluto 


En virtud de lo explicado en los párrafos an- 
teriores, la ecuación (31) y el principio en el 
cual está basada es importante desentrañar los 
mecanismos de formación de las plumas de 
trazador una vez que el soluto conservativo 
penetra en el agua. De esta forma, la concentra- 
ción pico inicial del soluto no es independiente 
de la concentración base en su entorno físico y 
tampoco es independiente del valor del caudal 
general. Es como si en el tiempo de formación 
de la pluma, su evolución estuviese predeter- 
minada por los parámetros generales del flujo, 
cuya información de alguna forma se transfiere 
a ese elemento de flujo; es muy probable que la 
naturaleza de esa transmisión sea de naturaleza 
electroquímica. La concentración de base Co es 
en principio “constante”, por lo tanto el ajuste 
requerido para obtener la concentración pico 
apropiada es función del tiempo; esto implica 
que su valor de ajuste debe ser calculado en un 
tiempo especial tc, como ya se dijo. 


La pendiente como dato clave de la condición 
de “equilibrio dinámico” 


Una observación cuidadosa de lo expresado en 
la introducción con respecto al papel prepon- 
derante de la pendiente del flujo como agente 
básico de las transformaciones energéticas en 
los cauces y por lo tanto en las condiciones de 
establecimiento del “equilibrio dinámico” que 
lleva al balance de masa de socavación y de 
deposición en el tramo indica que la función de 
estado q(t) no es independiente del mecanismo 
por el cual el flujo establece el valor preciso y 
específico de la pendiente, Sb, en el proceso de 
formación del cauce. 

Se debe decir que este proceso físico de 
fijación de la pendiente que actúa en específico 
para las condiciones de velocidad media, caudal 
y rugosidad tiene su símil en el proceso analítico 
para establecer cuál es el coeficiente longitudi- 
nal de dispersión del trazador (función de la 
pendiente), que vale realmente para el tramo 
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en su conjunto. En virtud de este símil se podría 
decir que la función de estado que sirve para la 
aplicación del principio de la equivalencia de 
las tasas de transporte es la misma que también 
se usa para la determinación de la pendiente. 

Para evidenciar esto es necesario hacer 
referencia a la ecuación de Elder y a su inter- 
pretación más reciente desde la teoría de los 
trazadores. 


Cálculo de la pendiente de un flujo en 
equilibrio dinámico según Elder 


En 1959, Elder propuso su ya célebre definición 
para el coeficiente longitudinal de dispersión, 
E (Elder, 1959); h es la profundidad media del 
flujo; Sb, la pendiente como fuerza generatriz 
gravitacional del movimiento del flujo, y g es la 
aceleración de la gravedad: 


E=5.93xhx,/hxgxS5, (32) 


Por supuesto, si la ecuación (21) se equipara 
con la ecuación (32), se debe aceptar que la 
definición de Elder de E debe ser también una 
función del tiempo. Es decir, puede moverse a 
través de la curva misma (figura 9). 


E(X|0 


Figura 9. E como función del tiempo. 
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Ubicación específica del coeficiente 
longitudinal de dispersión en la curva función 
del tiempo 


En realidad no basta afirmar que la definición de 
pendiente de Elder es una función del tiempo, 
sino también determinar en qué punto concreto 
de la curva estará situado el valor correspon- 
diente (punto negro en la curva de la figura 9). 
Para realizar esta tarea es necesario primero esti- 
mar el rango aproximado del tiempo a explorar 
(Constaín, 2014). 


Estimación del tiempo de acople “tc” 


Se parte de la ecuación (31), que da un pano- 
rama muy directo de las diferentes relaciones 
entre los parámetros de interés. Así, por ejem- 
plo, si se tiene en cuenta que q máximo puede 
valer 2.16, entonces la razón del numerador 
y del denominador tienen que ser del mismo 
orden. Además de esto, para ríos grandes en el 
trayecto inicial del trazador q >> q, que implica 
que C, >> Co. 


Cálculos de la pendiente, coeficiente de Chezy y 
función de estimación 


Una vez estimado el valor probable del tiempo 
que vale para aplicar la relación entre tasas de 
transporte, se procede a despejar la pendiente 
de la ecuación de Elder. Este procedimiento se 
debe hacer mediante “prueba y error”, para ir 
determinando para varios valores de E cuál es 
el que más se ajusta. El ajuste se hace mediante 
una función de estimación, F, que se desarrolla 
en seguida. Es prudente comenzar por aplicar 
el valor inicial encontrado mediante trazadores, 
según la ecuación (21). Si este primer intento no 
se ajusta, se usarán valores diferentes hasta lo- 
grar un ajuste apropiado, según se verá después 
en la aplicación experimental: 


E? 


$ (33) 
35.2xh x g 
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La función de estimación se establece igua- 
lando las dos definiciones de velocidad media 
del flujo, la nueva representada en la ecuación 
(20) y la clásica para flujo uniforme de Chezy, 
con R como radio hidráulico y Ch como el coefi- 
ciente de Chezy: 


A NI (34) 


QNYT 


A partir de allí se establece la función de 
estimación, F, que medirá la bondad de la 
aproximación que se haga para la pendiente 
mediante valores diversos de E. El ajuste es 
“ideal” cuando la función F es igual al coefi- 
ciente numérico de la ecuación de Elder: 


2 
Esq x0.215xt, (E 


x [2 5.93 (35) 
hxg 


Y, por lo tanto, en condición ideal: 


37.2 xh2 


A (36) 
É, xC?xRxa/S, 


Entonces, la función de estado q(t) que apa- 
rece en las ecuaciones (31) y (36) tiene el mismo 


valor que garantiza la validez del concepto de 
“equilibrio dinámico” en el tramo estudiado 
del cauce natural, ajustado de forma apropiada 
por la pendiente. La pendiente ajusta de modo 
simultáneo la dispersión y equivalencia de tasas 
de transporte. 

Cuando hay que realizar varias tentativas de 
“prueba y error” se debe recurrir de nuevo a la 
ecuación (21) para disponer de otros valores de 
E para determinar la pendiente con la debida 
aproximación. 


Aplicación del método a la 
determinación de la pendiente en un 
rio de montaña de mediano tamaño en 
Colombia 


Aspectos del cauce 


Se escogen los resultados experimentales con 
trazador salino realizados en el año 2007 en 
el río Pance en el Valle del Cauca, región sur 
occidental de Colombia. en En la figura 10 se 
muestran diversos aspectos del tramo medido. 


Instrumentación utilizada y aspectos del 
vertimiento de trazador iónico 


Para las mediciones con trazador iónico (sal co- 
mún) se utilizó un avanzado dispositivo deno- 
minado IDF (Inírida Deep Flow) desarrollado en 


Figura 10. Aspectos del cauce en el tramo medido. 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VIL, núm. 1, enero-febrero de 2016, pp. 89-105 


Tecnología y Ciencias del Agua, vol. VI, núm. 1, enero-febrero de 2016, pp. 89-105 


Constaín y Corredor, Universalidad de la función de estado que guía la dinámica de los solutos en los cauces naturales en... 


Colombia para operar en “tiempo real”, potente 
interfaz gráfica que permite un trabajo interac- 
tivo, así como guardar los aforos en memoria. 
En la figura 11 se muestran algunos aspectos 
de la instrumentación e inyección de trazador 
“aguas arriba”. En el experimento mostrado se 
usó una masa de sal común de 9 060 gramos 
disueltos inicialmente en agua, en un balde de 
unos 20 litros. 


Pantallas de la medición desde el equipo IDF 


En la figura 12 se tienen fotografías de las pan- 
tallas del instrumento con los datos respectivos 
de la medición. 

En la primera foto (izquierda superior) se 
muestra la curva experimental (trazo irregu- 
lar) y sobre puesto el modelo teórico logrado 
con la aplicación de la ecuación (33), dado los 
datos del cuadro 1. La cercanía del modelo 


teórico y la curva experimental es una garan- 
tía de la certeza de los datos obtenidos por el 
instrumento. 

La foto superior derecha corresponde al 
cálculo del caudal mediante el principio de 
conservación de la masa (área bajo la curva de 
trazador). La foto inferior izquierda correspon- 
de a la pantalla del instrumento en la que se 
muestra el cálculo de la función de estado qp(t) 
para la medición a X = 400 m. La foto inferior 
derecha muestra algunos datos calculados por 
el software del equipo. 


Cuadro de datos numéricos 


El cuadro 1 condensa las informaciones recogi- 
das por el IDF en la medición sobre el cauce. Se 
parte del dato de distancia del vertimiento y la 
medición que fue de X = 400 m para un ancho 
estimado de W = 20 metros. 


Figura 11. Aspectos del instrumento IDF y de la inyección de trazador. 
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Figura 12. Aspectos de las pantallas de las mediciones desde el equipo IDE. 


Desarrollo del estudio 
Estudio de tc mediante la relación de Elder 
Después de hacer varias tentativas en tiempos 


muy cercanos al instante del vertimiento, se 
mostrará en seguida que la dinámica correcta 


se puede caracterizar por el tiempo de acople tc 
= 2.6 s después del vertimiento, o sea a X=1.9 
m, si se considera que vale una velocidad media 
de U = 0.74 m/s. Los cálculos detallados para 
considerar que esto es correcto se muestran a 
continuación. 
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Cuadro 1. Resultados de trazador. 


RO 
E Masa, M Velocidad Caudal, | Concentración | Tiempo | Profundidad 
Vertimiento % 4d media, Ux DP Q pico Cp, t H 
e di m/s m/s Mgr/l Ss m 
Mgr/l 8 
X = 400 m 29.4 9.060 0.742 0.37 11.7 3.9 539 0.80 


Cálculo de Elder para la pendiente para tc =2.6 s 


Se inicia el proceso calculando el coeficiente 
longitudinal de dispersión, suponiendo que 
para este tiempo tan inmediato al vertimiento 


la función de estado q = 2.16: 


2 2 
O Pe xU ma 


A 2.16? x0.74?x0.215x2.6 
2 


E(t = 


=0.71m%/s (37) 


Se calcula luego la pendiente probable, supo- 
niendo que la profundidad media se conserva 
en todo el tramo considerado: 


E? 
O 
35.2x h” x9.81 


0.71 


 ———— ==00.0029 
35.2x0.8* x9.81 


(38) 


Se calcula luego el coeficiente aproximado 
de Chezy: 


u 0.74 


Cha a 
VRS  /0.8x0.0029 


=154m%/s (39) 


Se calcula luego la función de estimación, F: 


2 
PP xO 25M 5 


hxg 
15.4? 


a 0) 
2 N08x9.81 


Este valor se considera satisfactorio en tan- 
to es el ideal de 5.93 (coeficiente de Elder), y 


=2.167 x0.215x2.6x 
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entonces se pueden tomar como correctos los 
valores aproximados calculados antes. Se estima 
que este tiempo es el correcto para calcular la 
equivalencia de tasas de transporte. 


Cálculo de la razón de tasas de transporte atc=2.6 s 
Se empieza calculando de manera aproximada 
el coeficiente transversal de difusión, el ancho 


local de la pluma, para luego estimar el caudal 
local: 


e, =0.23xhx hxgx5 =0.23 x0.8 


x/0.8x 9.81 x 0.0029 =0.028 m/s (41) 
W, =,/3.22xe, xt =4/3.22x0.028x2.6 
= 0.48 m (42) 
Q, =W,xUxh=0.48x0.74x0.8 
=0.284 m/s = 284 1/s (43) 
La razón de caudales es entonces: 
9 _11700 4, 
q 284 q) 


Por lo tanto, la razón de concentraciones en 
estado estable a tc = 2.6 s es: 


q 
E. e 412 


= = =16.4 
C, 2.16x1.16 2.51 


(45) 


La concentración pico del trazador a tc = 2.6 
s después del vertimiento, con una masa de 9.06 
kg de soluto, es: 
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C, =16.4x 29.4 = 483 Mgr/1l (46) 


Verificación de los cálculos realizados mediante la 
determinación general de la curva de la Concentra- 
ción pico, Cp(t) 


Lo que se debe hacer ahora es construir esta fun- 
ción Cp(t) en función de los diferentes datos en 
instantes diferentes: A, en el punto de medición 
a t=539 s; B, en el punto de acople a tc=2.6 s. 

En el instante mismo de la inyección, at=0 
s. Para realizar la tarea completa se debe partir 
de lo explicado en el apartado “Modelo para la 
evolución de la concentración pico del trazador 
en el tiempo, Cp(t)”. Se comienza por conocer 


mMoÓS”, 


los distintos valores experimentales de “a”: 


A. Punto de medición at=539 s: 


0, mE 2 258 (47) 


É 
a, 2 29133 (48) 


m5” 


Como se dijo, uno puede constatar que “a. 
es —como q(t)— una función del tiempo. Por 
tal razón es necesario proponer un modelo ex- 
ponencial descendiente para este parámetro, de 
tal forma que el modelo coincida con los datos 
experimentales: 


PA 


C,(t) =913.3 09902113 (49) 


Se puede mostrar que para que la expresión 
ajuste, la constante “k” vale 0.00234, y entonces 
la concentración pico C,(t) reproduce de forma 
aproximada los datos experimentales para el 
tiempo de acople y medición (figura 13a y b). 


Verificación de la congruencia de los cálculos 


Tener el modelo general para C,(t) es muy 
ventajoso, pues se puede encontrar de forma 


aproximada la concentración pico del trazador 
en el momento mismo del vertimiento, el cual 
se considera tipo “slug” (súbito). 

De forma razonable se puede decir que la 
concentración pico inicial debe ser del mismo 
orden que la concentración de la mezcla inicial 
al momento de verter, es decir, la razón entre la 
masa de trazador vertido (9.060 kg) y el volu- 
men aproximado del balde de agua utilizado 
(20 litros): 


9 060 000 Mgr 


AA 
05) 201 


=453000 Mgr/1 (50) 

De acuerdo con la ecuación (49), se tiene que 
para un tiempo muy cercano a cero, la concen- 
tración pico del trazador es: 


C,( =0.0001 s) =913.3 002% + 


xE* = 424 651 Mgr/l (51) 


Es decir, un valor bastante cercano al espe- 
rado, con un error porcentual de 7%. Esta situa- 
ción se puede ver gráficamente en la figura 14. 


Conclusiones 


De lo analizado en este artículo se pueden 
identificar los siguientes aspectos sobre el ex- 
perimento de trazadores iónicos en el tramo del 
río Pance: 


1. Las aplicaciones experimentales de los linea- 
mientos teóricos presentados aquí fueron 
exitosas en la medida que los diferentes 
datos son altamente congruentes entre sí, 
aceptándose que representan —dentro del 
límite de las aproximaciones y simplifica- 
ciones usadas— una descripción fiable de la 
evolución del trazador y su relación con los 
parámetros geomorfológicos e hidráulicos. 


2. El concepto básico aquí discutido es el de 
“equivalencia de tasas de transporte” en 
los cauces naturales en la parte local y en 
la general; esto es una consecuencia de 
las propiedades estadísticas del flujo en 
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(b) Punto de medición 


Figura 13. Curva de Cp(t) replicando los datos experimentales. 
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Figura 13. Verificación aproximada para Cp(t = 0) en el punto de inyección. 
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condiciones estacionarias que aproximan 
el “equilibrio dinámico” analizado en de- 
talle por varios autores del siglo pasado, en 
especial Luna B. Leopold. Para un sistema 
fluvial cerrado, evolucionando de modo 
irreversible, no sólo se conserva la masa en 
sus fronteras sino que las tasas de transporte 
son iguales en cada punto del interior. Es un 
proceso en el cual la información global se 
incorpora en los fenómenos individuales. 


3. Esta equivalencia permite definir con 
precisión los mecanismos de formación 
de las plumas de trazador en los cauces 
naturales y por lo tanto de interés para los 
especialistas que estudian los procesos de 
contaminación en los estudios de calidad 
de aguas. El hecho de tener una informa- 
ción congruente en hidráulica, transporte 
dispersivo y geomorfología facilita mucho 
la calibración de los modelos a usarse de 
manera posterior. 
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